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Halo Shape and Relic Density Exclusions of Sommerfeld-Enhanced Dark Matter 

Explanations of Cosmic Ray Excesses 
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Dark matter with Sommerfeld-enhanced annihilation has been proposed to explain observed cos- 
mic ray positron excesses in the 10 GeV to TeV energy range. We show that the required enhance- 
ment implies thermal relic densities that are too small to be all of dark matter. We also show that the 
dark matter is sufficiently self-interacting that observations of elliptical galactic dark matter halos 
exclude large Sommerfeld enhancement for light force carriers. Resonant Sommerfeld enhancement 
does not modify these conclusions, and the astrophysical boosts required to resolve these discrep- 
ancies are disfavored, especially when significant self-interactions suppress halo substructure. 

PACS numbers: 95.35.+d, 95.85.Ry 



Introduction. Recently PAMELA [l(, ATIC (3, 
Fermi Q, and HESS Q have observed the spectrum of 
cosmic ray positrons with energies between 10 GeV and a 
few TeV. Some of the data show excesses over background 
expectations The excesses have plausible astrophys- 
ical explanations [1, 0]- At the same time, signals from 
many dark matter candidates are expected in this energy 
range, and this possibility has not escaped attention. 

By far the most researched possibility is that the ob- 
served positrons are produced by dark matter annihila- 
tion. If dark matter X is a thermal relic, the relic den- 
sity implies that its thermally-averaged annihilation cross 
section times relative velocity at freeze out is (<7 an i> rc i) ~ 
(Tq 11 = 3 x 10~ 26 cm 3 /s. Unfortunately, if this is the an- 
nihilation cross section now, the resulting signal is too 
small by two to three orders of magnitude to explain the 
observed cosmic ray excesses. 

A seemingly attractive solution is to postulate that 
dark matter interacts with a light force carrier <fi with fine 
structure constant ax = A 2 /(47r) For = 0, the 

annihilation is enhanced by the Sommerfeld factor [lfj 
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For massive </>, the enhancement is typically cutoff at a 
value oc arax/rn^ [1, 0, [H , 12 1 , but for fine-tuned choices 
of ax, rax and ra^,, there are also resonance regions 
where the enhancement may exceed this cutoff, as we dis- 
cuss below. The velocity of dark matter particles is ~ 1/3 
at freeze out and ~ 10~ 3 now. Sommerfeld enhancement 
therefore provides an elegant mechanism for boosting an- 
nihilations now. Constraints from dark matter annihila- 
tion in protohalos with v Te \ ~ 10 -8 exclude = 
However, taking mx ~ TeV and ~ MeV — GeV, and 
assuming (er an i> rG i) s» (Tg h , one may still generate S ~ 10 3 
to explain the positron excesses, while the cutoff allows 
one to satisfy the protohalo constraint. 

Of course, for a viable solution, dark matter must not 
only annihilate with the correct rate, it must also be pro- 
duced with the right density and form structure in accord 



with observations. Here we find that the desired thermal 
relic density cannot be achieved in Sommerfeld-enhanced 
models designed to explain the positron excesses. In ad- 
dition, we show that the new force carrier </> induces dark 
matter self-interactions that may contradict current ob- 
servations. As is well-known, if <p were massless, the 
resulting long range Coulomb force would lead to large 
energy transfers that make halos spherical, and observa- 
tions of triaxial halos constrain this possibility 
For to^ ~ 100 MeV, the force's range is only ~ 10 fm, 
but, as we show below, the implied cross section is still 
large enough to play a role in galactic dynamics. 

Thermal Relic Density. If XX annihilation is en- 
hanced by <p exchange, there is an "irreducible" anni- 
hilation process A" A" — > <j>4> through ^-channel X. For 
rrirf, -C rax, the thermally-averaged annihilation cross 
section is the typical WIMP cross section 
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with 0(1) corrections depending on the details of the 
initial and final states. Requiring that (cr^v^i) be small 
enough that X can be all of the dark matter implies 



ax < ^/cr^ h /7r ra x 



(3) 



This bound is conservative. In fact, the Sommerfeld ef- 
fect enhances the annihilation cross section even at freeze 
out d, [lU , and the bound may be significantly strength- 
ened in the presence of other annihilation channels. 

Self-Interactions. Self-interactions allow dark matter 
particles to transfer energy. The average rate for dark 
matter particles to change velocities by 0(1) factors is 

r fc = J d 3 v 1 d 3 v 2 f(v 1 )f(v 2 ) (nxv^iar) (v?Jv$) , (4) 

where f(v) — e~ y2 ^ v o /(vo^/n) 3 is the dark matter's as- 
sumed (Maxwellian) velocity distribution, nx is its num- 
ber density, Vroi = \vi— v 2 \, and ot = J dCl* (da/dn*)(l — 
cos 9*) is the energy transfer cross section, where 0* is the 
scattering angle in the center-of-mass frame. 
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Dark matter particles coupled to a massive force car- 
rier scatter through the Yukawa potential V(r) = 
—axe~ m4,r /r. In the Born approximation, keeping only 
the dominant i-channel contribution present in all inter- 
actions, the transfer cross section is 
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where (3 = 2ax'm c f > / (mxv 2 el ) is the ratio of the potential 
energy at r ~ 1 to the kinetic energy of the particle, 
and R = mxv^/m^ is the ratio of the interaction range 
to the dark matter particle's de Broglie wavelength. For 
typical values of interest here, i> rc i ~ 10 -3 and rax /m<f, > 

10 3 , and so R > 1. For E> 1, a T « J?^- (lni? 2 - l). 

As in the Coulomb case, this is greatly enhanced for small 
v re i, but here the finite interaction length of the Yukawa 
potential cuts off the logarithmic divergence. 

Equation ([5]) receives significant corrections in the 
strong interaction regime, where /3 » 1. Our focus in 
this work will be on the R^> 1 region of parameter space. 
In this region, quantum effects are subdominant and so 
classical studies of particles moving in Yukawa po- 
tentials are applicable. Although the authors of these 
studies were interested in slow and highly charged par- 
ticles moving in plasmas with screened Coulomb poten- 
tials, they approximated these potentials by Yukawa po- 
tentials, and so their results are exactly applicable in the 
current context. The numerical results of these studies 
are accurately reproduced by 
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, 0.1 < $ < 1000 . (6) 



We use these analytical fits to obtain the results below. 

Halo Shapes. Self-interactions that are strong enough 
to create 0(1) changes in the energies of dark matter 
particles will isotropize the velocity dispersion and cre- 
ate spherical halos. These expectations are borne out by 
simulations of self-interacting dark matter in the hard 
sphere limit 18|-2(|. The shapes of dark matter halos 



of elliptical galaxies and clusters are decidedly elliptical, 
which constrains self-interactions [2l|. The ellipticity of 
galactic halos provides the strongest constraints on these 
models [15j . To implement these constraints, we consider 
the well-studied, nearby (about 25 Mpc away) elliptical 
galaxy NGC 720. In Ref. (22j |. X-ray isophotes were used 
to extract the ellipticity of the underlying matter dis- 
tribution. Comparing it to the ellipticity induced by the 
stellar mass profile, the dark matter halo of NGC 720 was 
found to be elliptical at about 5 kpc and larger radii. 

To compute Tk, we use the measured total mass pro- 
file and the decomposition into stars plus dark matter for 
NGC 720 23] and obtain the radial velocity dispersion 



v 2 (r) — «g(r)/2 and the dark matter density. For the 
radius we pick 5 kpc. Our constraints would be stronger 
if we could use the higher densities inside this radius, but 
the constraints on the ellipticity weaken for radii below 
5 kpc [22^ . For the dark matter density, we choose the 
average value within 5 kpc, which is roughly 4 GeV/cm 3 . 
To compute the dispersion, we assume isotropy and that 
the total (stellar plus dark matter) mass profile scales 
approximately linearly with radius. For an NFW profile 
with best fit scale radius (23|, v 2 (r) ~ (240 km/s) 2 . Vary- 
ing within the quoted error range for the scale radius 23] 
only changes this dispersion by about 10%. 

Results. To derive constraints on the particle physics 
parameters from the observed halo shapes, we require 



r; 



1 > 10 10 years 



(7) 



i.e., that the average time for self- interactions to cre- 
ate 0(1) changes in dark matter particle velocities is 
greater than the galaxy's lifetime. Imposing Eqs. ([3]) 
and ((7|) from the relic density and the observation of el- 
lipticity in the dark matter halo of NGC 720 yields the 
constraints shown in Fig. [TJ The relic density constraint 
is independent of and the extremely stringent halo 
shape constraint for = [l5| remains significant for 
m 4> ^5 30 MeV. The crucial point is that when the inter- 
action range is larger than the de Broglie wavelength, al- 
though the Coulomb logarithm enhancement is lost, the 
enhancement from low u re i remains. Note that our as- 
sumption of a locally Gaussian velocity distribution is 
supported by recent simulation of Milky Way-sized dark 
matter halos [2§] . does not change by more than a fac- 
tor of about 2 when we allow the distribution to become 
anisotropic or introduce a velocity cut-off at the escape 
speed. At the same time, we have checked that our halo 
shape bounds are consistent with the predictions from 
simulations with hard sphere scattering fljj . 

In Fig. [2] we present the regions of the (mx , S) plane 
required to explain PAMELA and Fermi as determined 
241 . These are for = 250 MeV, which is 



Ref. 



large enough to allow contributions to positrons through 
<p — > fi + fi~ , but small enough to forbid contributions to 
anti-protons, where no excess is seen 0. Upper bounds 
from relic density and halo shapes are also given. We 
see that the large Sommerfeld enhancements required 
to explain the positron excesses are significantly ex- 
cluded by the relic density constraint for all mx- For 
m.$ ^ 30 MeV, the halo shape constraints also exclude 
the required Sommerfeld enhancements. 

Discussion. The results of Fig. [2] are not surpris- 
ing. For the relic density, the WIMP miracle implies 
that for mx ~ 250 GeV, the correct relic density is 
obtained for a ~ 10~ 2 . Given u re i ~ 10~ 3 , this im- 
plies an upper bound of S ~ 10, and this bound scales 
as mx- Of course, X need not be all the dark matter, 
but in this case, the Sommerfeld-enhanced flux scales as 
n 2 (o~ an v)S ~ a x l , and so the signal is maximal for S ~ 1. 
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FIG. 1: Regions above the contours are excluded by the relic 
density constraint and by halo ellipticity observations for the 
indicated. The classical approximation used to obtain the 
halo bounds becomes inaccurate for m$ > 100 MeV. 
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FIG. 2: Upper bounds on Sommerfeld enhancement fac- 
tor S from relic density (solid), along with PAMELA- and 
Fermi- favored re gion s and the best fit point (nix,S) = 
(2.35 TeV, 1500) [241, all for m # = 250 MeV. Halo shape 
bounds are also shown for the values of indicated (dashed). 



In deriving our results, we have ignored the cutoff of 
the Sommerfeld enhancement factor for massive 4>. In- 
cluding this cutoff will reduce the maximal possible S for 
low mx , strengthening the disagreement between the al- 
lowed values of S and the experimentally favored regions. 
To reduce the disagreement, one might consider resonant 
Sommerfeld enhancement. As with resonances from addi- 



tional postulated particles [25[, these resonances require 
fine-tuning and are bounded by astrophysical observa- 
tions 26]. In addition, resonance enhancement occurs 
at m^/mx — Gax/i^n 2 ), n = 1,2,... [27j and is sig- 
nificant only for low n. For m^ ~ GeV and the rele- 
vant range of ax <L 0.01, this implies mx ^ 500 GeV; 
the resonances are, therefore, ineffective in reaching the 



favored regions given in Fig. [2j Most importantly, as 
noted above, our bounds are conservative in that they 
do not include the Sommerfeld effect on freeze out [16j . 
This effect suppresses the largest possible S, especially at 
resonances. Self-consistcntly including the effects of res- 
onances on annihilation in both the early Universe and 
now, we find that the maximal possible enhancement fac- 
tor is S ~ 100, even allowing for resonances [28j . 

As an alternative approach to evade the relic density 
constraints, one may consider other production mecha- 
nisms or modify early Universe cosmology, but this sac- 
rifices the WIMP miracle and also removes the moti- 
vation for considering Sommerfeld enhancement in the 
first place. Alternatively, one might appeal to boosts of 
~ 10 from cold and dense dark matter substructure in 
the local neighborhood. Such large values at a distance 
of only 10 kpc from the Milky Way center are, however, 
not motivated by simulations with collisionless dark mat- 
ter [29[ . The presence of the stellar disk would further re- 
duce these expectations. The self-scatterings among the 
particles in the substructure would also serve to reduce 
the inner densities [HI and hence the expected boost. In 
addition, for m^ < 30 MeV, interactions with the dark 
matter particles of the Milky Way could evaporate sub- 
structure because i> r ei is much larger than the internal 
velocity dispersion of the substructure [3(| . 

The halo shape bounds are obtained from inferred dark 
matter halo ellipticity, which depends on merger histories 
and the environment. For example, a major merger at a 
redshift of z = 0.5 for NGC 720 would effectively halve 
the age that TjT 1 should be compared to and weaken the 
bound on m^ by roughly a factor of \[2. However, the 
lack of large scale disturbances in the gas argues against 
such a recent major merger. These bounds may be made 
more robust by deeper data sets of NGC 720, which will 
further constrain point source contamination and rota- 
tion or large scale disturbances in the g well as by 
measuring ellipticities and mass profiles in other galaxies 
and clusters [32j. 

A second prediction of strongly self-interacting dark 
matter is the formation of constant density cores, if 
gravo-thermal collapse does not occur. The time scale 
for the formation of these cores is of order F^ , suggest- 
ing that NGC 720 should have O(kpc) sized core. Fu- 
ture tests for the presence of cores in galaxy and cluster 
halos may provide comparable or stronger limits. Self- 
interactions should also dramatically alter the dark mat- 
ter halos of smaller galaxies, such as the dwarf galax- 
ies in the Local Group. The central dark matter densi- 
ties measured in these dwarf satellites of the Milky Way 
are C(GeV/cm 3 ) [33| and fit neatly within the standard 
CDM predictions. For the parameter space disfavore d by 
NGC 720 observations, and using simulation results [18| . 
we estimate that core sizes would be of order the lumi- 
nous extent of the dwarfs or larger. The tidal force of the 
Milky Way would significantly reduce the central densi- 
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ties of the dwarfs with such large cores and likely make 
it impossible to explain the large observed dark matter 
densities in all the dwarfs 3J) . In parameter regions with 
more moderate Sommerfeld enhancements, these cores 
would be smaller and consistent with current data 311 ] . 

Conclusions. Cosmic positron data have motivated 
dark matter candidates with Sommerfeld-enhanced an- 
nihilations. The required enhancement is large, requir- 
ing large couplings to light force carriers. Annihilation 
to these force carriers provides an upper limit on the 
thermal relic abundance of these dark matter candidates. 
With or without resonances, this constraint excludes the 
existence of enhancements that can explain the positron 
excesses. These models also predict self-interactions that 
may make galactic dark matter halos spherical. The el- 
lipticity of the halo of NGC 720 also excludes the required 
Sommerfeld enhancements for < 30 MeV. Interest- 
ingly, viable models with moderate Sommerfeld enhance- 
ments, although unable to explain the positron data, may 
predict constant density spherical cores in small galactic 
halos and other departures from the standard cold dark 
matter paradigm that are consistent with current data. 
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Note added. As the first version of this work was being 
completed, we learned of related work in progress. This 
work [3^1 agrees with Eq. (j6|) in the classical regime. 
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